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PERCEPTUAL CODING OF AUDIO SIGNALS USING SEPARATED 
IRRELEVANCY REDUCTION AND REDUNDANCY REDUCTION 

Cross-Reference to Related Applications 

The present invention is related to United States Patent Application entitled 
"Method and Apparatus for Representing Masked Thresholds in a Perceptual Audio Coder," 
(Attorney Docket Number Edler 2-2-6), United States Patent Application entitled "Perceptual 
Coding of Audio Signals Using Cascaded Filterbanks for Performing Irrelevancy Reduction and 
Redundancy Reduction With Different Spectral/Temporal Resolution," (Attorney Docket 
Number Edler 3-4), United States Patent Application entitled "Method and Apparatus for 
Reducing Aliasing in Cascaded Filter Banks," (Attorney Docket Number Schuller 5) and United 
States Patent Application entitled "Method and Apparatus for Detecting Noise-Like Signal 
Components," (Attorney Docket Number Fink Faller 3), filed contemporaneously herewith, 
assigned to the assignee of the present invention and incorporated by reference herein. 



Field of the Invention 

The present invention relates generally to audio coding techniques, and more 
particularly, to perceptually-based coding of audio signals, such as speech and music signals. 



Background of the Invention 

Perceptual audio coders (PAC) attempt to minimize the bit rate requirements for 
the storage or transmission (or both) of digital audio data by the application of sophisticated 
hearing models and signal processing techniques. Perceptual audio coders (PAC) are described, 
for example, in D. Sinha et al., "The Perceptual Audio Coder," Digital Audio, Section 42, 42-1 to 
42-18, (CRC Press, 1998), incorporated by reference herein. In the absence of channel errors, a 
PAC is able to achieve near stereo compact disk (CD) audio quality at a rate of approximately 
128 kbps. At a lower rate of 96 kbps, the resulting quality is still fairly close to that of CD audio 
for many important types of audio material. 

Perceptual audio coders reduce the amount of information needed to represent an 
audio signal by exploiting human perception and minimizing the perceived distortion for a given 




Edler 1-4 



bit rate. Perceptual audio coders first apply a time-frequency transform, which provides a 
compact representation, followed by quantization of the spectral coefficients. FIG. 1 is a 
schematic block diagram of a conventional perceptual audio coder 100. As shown in FIG. 1, a 
typical perceptual audio coder 100 includes an analysis filterbank 110, a perceptual model 120, a 
5 quantization and coding block 130 and a bitstream encoder/multiplexer 140. 

The analysis filterbank 110 converts the input samples into a sub-sampled spectral 
representation. The perceptual model 120 estimates the masked threshold of the signal. For each 
spectral coefficient, the masked threshold gives the maximum coding error that can be introduced 
into the audio signal while still maintaining perceptually transparent signal quality. The 
10 quantization and coding block 130 quantizes and codes the prefilter output samples according to 
the precision corresponding to the masked threshold estimate. Thus, the quantization noise is 
^ hidden by the respective transmitted signal. Finally, the coded prefilter output samples and 
yo additional side information are packed into a bitstream and transmitted to the decoder by the 
gj bitstream encoder/multiplexer 140. 

jjfj; 15 FIG. 2 is a schematic block diagram of a conventional perceptual audio decoder 

Si 200. As shown in FIG. 2, the perceptual audio decoder 200 includes a bitstream 
3 ' decoder/demultiplexer 210, a decoding and inverse quantization block 220 and a synthesis 
y filterbank 230. The bitstream decoder/demultiplexer 210 parses and decodes the bitstream 
O yielding the coded prefilter output samples and the side information. The decoding and inverse 
p 20 quantization block 220 performs the decoding and inverse quantization of the quantized prefilter 

output samples. The synthesis filterbank 230 transforms the prefilter output samples back into 

the time-domain. 

Generally, the amount of information needed to represent an audio signal is 
reduced using two well-known techniques, namely, irrelevancy reduction and redundancy 
25 removal. Irrelevancy reduction techniques attempt to remove those portions of the audio signal 
that would be, when decoded, perceptually irrelevant to a listener. This general concept is 
described, for example, in U.S. Pat. No. 5,341,457, entitled "Perceptual Coding of Audio 
Signals," by J. L. Hall and J. D. Johnston, issued on Aug. 23, 1994, incorporated by reference 
herein. 
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Currently, most audio transform coding schemes implemented by the analysis 
filterbank 110 to convert the input samples into a sub-sampled spectral representation employ a 
single spectral decomposition for both irrelevancy reduction and redundancy reduction. The 
redundancy reduction is obtained by dynamically controlling the quantizers in the quantization 
5 and coding block 130 for the individual spectral components according to perceptual criteria 
contained in the psychoacoustic model 120. This results in a temporally and spectrally shaped 
quantization error after the inverse transform at the receiver 200. As shown in FIGS. 1 and 2, the 
psychoacoustic model 120 controls the quantizers 130 for the spectral components and the 
corresponding dequantizer 220 in the decoder 200. Thus, the dynamic quantizer control 

10 information needs to be transmitted by the perceptual audio coder 100 as part of the side 
information, in addition to the quantized spectral components. 

The redundancy reduction is based on the decorrelating property of the transform. 
For audio signals with high temporal correlations, this property leads to a concentration of the 
signal energy in a relatively low number of spectral components, thereby reducing the amount of 

15 information to be transmitted. By applying appropriate coding techniques, such as adaptive 
Huffman coding, this leads to a very efficient signal representation. 

One problem encountered in audio transform coding schemes is the selection of 
the optimum transform length. The optimum transform length is directly related to the frequency 
resolution. For relatively stationary signals, a long transform with a high frequency resolution is 

20 desirable, thereby allowing for accurate shaping of the quantization error spectrum and providing 
a high redundancy reduction. For transients in the audio signal, however, a shorter transform has 
advantages due to its higher temporal resolution. This is mainly necessary to avoid temporal 
spreading of quantization errors that may lead to echoes in the decoded signal. 

As shown in FIG. 1, however, conventional perceptual audio coders 100 typically 

25 use a single spectral decomposition for both irrelevancy reduction and redundancy reduction. 
Thus, the spectral/temporal resolution for the redundancy reduction and irrelevancy reduction 
must be the same. While high spectral resolution yields a high degree of redundancy reduction, 
the resulting long transform window size causes reverbation artifacts, impairing the irrelevancy 
reduction. A need therefore exists for methods and apparatus for encoding audio signals that 

30 permit independent selection of spectral and temporal resolutions for the redundancy reduction 
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and irrelevancy reduction. A further need exists for methods and apparatus for encoding speech 
as well as music signals using a psychoacoustic model (a noise-shaping filter) and a transform. 



Summary of the Invention 

5 Generally, a perceptual audio coder is disclosed for encoding audio signals, such 

as speech or music, with different spectral and temporal resolutions for the redundancy reduction 
and irrelevancy reduction. The disclosed perceptual audio coder separates the psychoacoustic 
model (irrelevancy reduction) from the redundancy reduction, to the extent possible. The audio 
signal is initially spectrally shaped using a prefilter controlled by a psychoacoustic model. The 
10 prefilter output samples are thereafter quantized and coded to minimize the mean square error 
(MSE) across the spectrum. 

According to one aspect of the invention, the disclosed perceptual audio coder 

3 

0 uses fixed quantizer step-sizes, since spectral shaping is performed by the pre-filter prior to 
& quantization and coding. Thus, additional quantizer control information does not need to be 

L; 15 transmitted to the decoder, thereby conserving transmitted bits. 

3 

4 The disclosed pre-filter and corresponding post-filter in the perceptual audio 

decoder support the appropriate frequency dependent temporal and spectral resolution for 
!f irrelevancy reduction. A filter structure based on a frequency-warping technique is used that 
3 allows filter design based on a non-linear frequency scale. 

3 20 The characteristics of the pre-filter may be adapted to the masked thresholds (as 

^ generated by the psychoacoustic model), using techniques known from speech coding, where 
linear-predictive coefficient (LPC) filter parameters are used to model the spectral envelope of 
the speech signal. Likewise, the filter coefficients may be efficiently transmitted to the decoder 
for use by the post-filter using well-established techniques from speech coding, such as an LSP 
25 (line spectral pairs) representation, temporal interpolation, or vector quantization. 

A more complete understanding of the present invention, as well as further 
features and advantages of the present invention, will be obtained by reference to the following 
detailed description and drawings. 



-4- 




Edler 1-4 



Brief Description of the Drawings 

FIG. 1 is a schematic block diagram of a conventional perceptual audio coder; 

FIG. 2 is a schematic block diagram of a conventional perceptual audio decoder 
corresponding to the perceptual audio coder of FIG. 1; 

FIG. 3 is a schematic block diagram of a perceptual audio coder according to the 
present invention and its corresponding perceptual audio decoder; 

FIG. 4. illustrates an FIR predictor of order P, and the corresponding IIR 

predictor; 

FIG. 5 illustrates a first order allpass filter; and 

FIG. 6 is a schematic diagram of an FIR filter and a corresponding IIR filter 
exhibiting frequency warping in accordance with one embodiment of the present invention. 

Detailed Description 

FIG. 3 is a schematic block diagram of a perceptual audio coder 300 according to 
the present invention and its corresponding perceptual audio decoder 350, for communicating an 
audio signal, such as speech or music. While the present invention is illustrated using audio 
signals, it is noted that the present invention can be applied to the coding of other signals, such as 
the temporal, spectral, and spatial sensitivity of the human visual system, as would be apparent to 
a person of ordinary skill in the art, based on the disclosure herein. 

According to one feature of the present invention, the perceptual audio coder 300 
separates the psychoacoustic model (irrelevancy reduction) from the redundancy reduction, to the 
extent possible. Thus, the perceptual audio coder 300 initially performs a spectral shaping of the 
audio signal using a prefilter 310 controlled by a psychoacoustic model 315. For a detailed 
discussion of suitable psychoacoustic models, see, for example, D. Sinha et al, "The Perceptual 
Audio Coder," Digital Audio, Section 42, 42-1 to 42-18, (CRC Press, 1998), incorporated by 
reference above. Likewise, in the perceptual audio decoder 350, a post-filter 380 controlled by 
the psychoacoustic model 315 inverts the effect of the pre-filter 310. As shown in FIG. 3, the 
filter control information needs to be transmitted in the side information, in addition to the 
quantized samples. 
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Quantizer/Coder 

The prefilter output samples are quantized and coded at stage 320. As discussed 
further below, the redundancy reduction performed by the quantizer/coder 320 minimizes the 
mean square error (MSE) across the spectrum. 

Since the pre-filter 310 performs spectral shaping prior to quantization and 
coding, the quantizer/coder 320 can employ fixed quantizer step-sizes. Thus, additional 
quantizer control information, such as individual scale factors for different regions of the 
spectrum, does need not need to be transmitted to the perceptual audio decoder 350. 

Well-known coding techniques, such as adaptive Huffman coding, may be 
employed by the quantizer/coder stage 320. If a transform coding scheme is applied to the pre- 
filtered signal by the quantizer/coder 320, the spectral and temporal resolution can be fully 
optimized for achieving a maximum coding gain under a mean square -error (MSE) criteria. As 
discussed below, the perceptual noise shaping is performed by the post-filter 380. Assuming the 
distortions introduced by the quantization are additive white noise, the temporal and spectral 
structure of the noise at the output of the decoder 350 is fully determined by the characteristics of 
the post-filter 380. It is noted that the quantizer/coder stage 320 can include a filterbank such as 
the analysis filterbank 110 shown in FIG. 1. Likewise, the decoder/dequantizer stage 360 can 
include a filterbank such as the synthesis filterbank 230 shown in FIG. 2. 

Pre-Filter/Post-Filter Based on Psvchoacoustic Model 

One implementation of the pre-filter 310 and post-filter 380 is discussed further 
below in a section entitled "Structure of the Pre-Filter and Post-Filter." As discussed below, it is 
advantageous if the structure of the pre-filter 310 and post-filter 380 also supports the appropriate 
frequency dependent temporal and spectral resolution. Therefore, a filter structure based on a 
frequency-warping technique is used which allows filter design on a non-linear frequency scale. 

For using the frequency warping technique, the masked threshold needs to be 
transformed to an appropriate non-linear (i.e. warped) frequency scale as follows. Generally, the 
resulting procedure to obtain the filter coefficients g is: 

- Application of the psychoacoustic model gives a masked threshold as power 
(density) over frequency. 
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- A non-linear transformation of the frequency scale according to the frequency 
warping, as discussed below, gives a transformed masked threshold. 

- Application of LPC analysis / modeling techniques leads to LPC filter 
coefficients h, which can be quantized and coded using a transformation to lattice coefficients or 

5 LSPs 

- for use in the warped filter structure shown in FIG. 6, the LPC filter coefficients, 
h, need to be converted to filter coefficients, g 

The characteristics of the filter 310 may be adapted to the masked thresholds (as 
generated by the psychoacoustic model 315), using techniques known from speech coding, where 
10 linear-predictive coefficient (LPC) filter parameters are used to model the spectral envelope of 
the speech signal. In conventional speech coding techniques, the LPC filter parameters are 
usually generated in a way that the spectral envelope of the analysis filter output signal is 

S maximally flat. In other words, the magnitude response of the LPC analysis filter is an 

I n 

approximation of the inverse of the input spectral envelope. The original envelope of the input 
5* 15 spectrum is reconstructed in the decoder by the LPC synthesis filter. Therefore, its magnitude 
Sj response has to be an approximation of the input spectral envelope. For a more detailed 
5 " discussion of such conventional speech coding techniques, see, for example, W.B. Kleijn and 
5 K.K. Paliwal, "An Introduction to Speech Coding," in Speech Coding and Synthesis, 
Q Amsterdam: Elsevier (1995), incorporated by reference herein. 

^ 20 Similarly, the magnitude responses of the psychoacoustic post-filter 380 and pre- 

B filter 310 should correspond to the masked threshold and its inverse, respectively. Due to this 
similarity, known LPC analysis techniques can be applied, as modified herein. Specifically, the 
known LPC analysis techniques are modified such that the masked thresholds are used instead of 
short-term spectra. In addition, for the pre-filter 310 and the post-filter 380, not only the shape of 
25 the spectral envelope has to be addressed, but the average level has to be included in the model as 
well. This can be achieved by a gain factor in the post-filter 380 that represents the average 
masked threshold level, and its inverse in the pre-filter 310. 

Likewise, the filter coefficients may be efficiently transmitted using well- 
established techniques from speech coding, such as an LSP (line spectral pairs) representation, 
30 temporal interpolation, or vector quantization. For a detailed discussion of such speech coding 
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techniques, see, for example, F.K. Soong and B.-H. Juang, "Line Spectrum Pair (LSP) and 
Speech Data Compression," in Proc. ICASSP (1984), incorporated by reference herein. 

One important advantage of the pre-filter concept of the present invention over 
standard transform audio coding techniques is the greater flexibility in the temporal and spectral 
5 adaptation to the shape of the masked threshold. Therefore, the properties of the human auditory 
system should be taken into account in the selection of the filter structures. For a more detailed 
discussion of the characteristics of the masking effects, see, for example, M. R. Schroeder et aL, 
"Optimizing Digital Speech Coders By Exploiting Masking Properties Of The Human Ear," 
Journal of the Acoust. Soc. Am., v. 66, 1647-1652 (Dec. 1979); and J. H. Hall, "Auditory 
10 Psychophysics For Coding Applications," The Digital Signal Processing Handbook (V. Madisetti 
and D. B. Williams, eds.), 39-1:39-22, CRC Press, IEEE Press (1998), each incorporated by 
_ reference herein. 

j] Generally, the temporal behavior is characterized by a relatively short rise time 

^ even starting before the onset of a masking tone (masker) and a longer decay after it is switched 

ff 1 15 off. The actual extent of the masking effect also depends on the masker frequency leading to an 
\| increase of the temporal resolution with increasing frequency. 

\* For stationary single tone maskers, the spectral shape of the masked threshold is 

2 spread around the masker frequency with a larger extent towards higher frequencies than towards 

Q lower frequencies. Both of these slopes strongly depend on the masker frequency leading to a 

fl, 20 decrease of the frequency resolution with increasing masker frequency. However, on the non- 
El linear "Bark scale," the shapes of the masked thresholds are almost frequency independent. This 
Bark scale covers the frequency range from zero (0) to 20 kHz with 24 units (Bark). 

While these characteristics have to be approximated by the psychoacoustic model 
315, it is advantageous if the structure of the pre- filter 310 and post- filter 380 also supports the 
25 appropriate frequency dependent temporal and spectral resolution. Therefore, as previously 
indicated, the selected filter structure described below is based on a frequency-warping technique 
that allows filter design on a non-linear frequency scale. 

Structure of the Pre-Filter and Post-Filter 
The pre-filter 310 and post- filter 380 must model the shape of the masked 
30 threshold in the decoder 350 and its inverse in the encoder 300. The most common forms of 
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predictors use a minimum phase finite-impulse response (FIR) filter in the encoder 300 leading 
to an IIR filter in the decoder. FIG. 4. illustrates an FIR predictor 400 of order P, and the 
corresponding IIR predictor 450. The structure shown in FIG. 4 can be made time-varying quite 
easily, since the actual coefficients in both filters are equal and therefore can be modified 
synchronously. 

For modeling masked thresholds, a representation with the capability to give more 
detail at lower frequencies is desirable. For achieving such an unequal resolution over frequency, 
a frequency-warping technique, described, for example, in H. W. Strube, "Linear Prediction on a 
Warped Frequency Scale," J. of the Acoust. Soc. Am., vol. 68, 1071-1076 (1980), incorporated 
by reference herein, can be applied effectively. This technique is very efficient in the sense of 
achievable approximation accuracy for a given filter order which is closely related to the required 
amount of side information for adaptation. 

Generally, the frequency-warping technique is based on a principle which is 
known in filter design from techniques like lowpass-lowpass transform and lowpass-bandpass 
transform. In a discrete time system an equivalent transformation can be implemented by 
replacing every delay unit by an all-pass. A frequency scale reflecting the non-linearity of the 
"critical band" scale would be the most appropriate. See, M. R. Schroeder et al., "Optimizing 
Digital Speech Coders By Exploiting Masking Properties Of The Human Ear," Journal of the 
Acoust. Soc. Am., v. 66, 1647-1652 (Dec. 1979); and U. K. Laine et al., "Warped Linear 
Prediction (WLP) in Speech and Audio Processing," in IEEE Int. Conf Acoustics, Speech, 
Signal Processing, EH-349 - 1H-352 (1994), each incorporated by reference herein. 

Generally, the use of a first order allpass filter 500, shown in FIG. 5, gives a 
sufficient approximation accuracy. However, the direct substitution of the first order allpass 
filter 500 into the FIR 400 of FIG. 4 is only possible for the pre-filter 310. Since the first order 
allpass filter 500 has a direct path without delay from its input to the output, the substitution of 
the first order allpass filter 500 into the feedback structure of the IIR 450 in FIG. 4 would result 
in a zero-lag loop. Therefore, a modification of the filter structure is required. In order to allow 
synchronous adaptation of the filter coefficients in the encoder and decoder, both systems should 
be modified as described hereinafter. 
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In order to overcome this zero-lag problem, the delay units of the original 
structure (FIG. 4) are replaced by first order IIR filters containing only the feedback part of the 
first order allpass filter 500, as described in H.W. Strube, incorporated by reference above. FIG. 
6 is a schematic diagram of an FIR filter 600 and an IIR filter 650 exhibiting frequency warping 
in accordance with one embodiment of the present invention. The coefficients of the filter 600 
need to be modified to obtain the same frequency as a structure with allpass units. The 
coefficients, g k (0 <k <P), are obtained from the original LPC filter coefficients with the 
following transformation: 



8t =jlC h ,h n withC hl =f"l(l-a 2 )*(-ar 

n=k W 



10 The use of a first order allpass in the FIR filter 600 leads to the following mapping of the 
frequency scale: 

asmco 

xd - co + arctan 

l-acos<y 

The derivative of this function: 

dm \-a 2 



v(a>) = - 



dco \ + a 2 -2acosco 

15 indicates whether the frequency response of the resulting filter 600 appears compressed (v > 1) or 
stretched (v < 1). The warping coefficient a should be selected depending on the sampling 
frequency. For example, at 32 kHz, a warping coefficient value around 0.5 is a good choice for 
the pre-filter application. 

It is noted that the pre-filter method of the present invention is also useful for 

20 audio file storage applications. In an audio file storage application, the output signal of the pre- 
filter 310 can be directly quantized using a fixed quantizer and the resulting integer values can be 
encoded using lossless coding techniques. These can consist of standard file compression 
techniques or techniques highly optimized for lossless coding of audio signals. This approach 
opens the applicability of techniques that, up to now, were only suitable for lossless compression 

25 towards perceptual audio coding. 

It is to be understood that the embodiments and variations shown and described 
herein are merely illustrative of the principles of this invention and that various modifications 
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may be implemented by those skilled in the art without departing from the scope and spirit of the 
invention. 



